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Abstract. 

Standard solar models predict a solar luminosity that gradually increased by about 
30% over the past 4.5 billion years. Under the faint sun, Earth should have been frozen 
solid for most of its existence. Yet, running water is observed to have been present since 
very early in Earth's history. This enigma is known as the faint sun paradox. We show 
here that it can be partially resolved once we consider the cooling effect that cosmic rays 
are suspected to have on the global climate and that the younger sun must have had 
a stronger solar wind, such that it was more effective at stopping cosmic rays from reach- 
ing Earth. The paradox can then be completely resolved with the further contribution 
of modest greenhouse gas warming. When we add the cosmic ray flux modulation by a 
variable star formation rate in the Milky Way, we recover the long term glacial activ- 
ity on Earth. As to the future, we find that the average global temperature will increase 
by typically 10°K in the coming 2 Gyr. 



1. Introduction 

According to standard solar models, the solar luminos- 
ity increased from about 70% of the present solar luminos- 
ity at 4.5 Gyr before present (BP) to its present value. If 
Earth were a black body, its temperature would have been 
~ 25°K lower, enough to have kept large parts of it frozen 
until about 1-2 Gyr BP. Besides however the past Eon, and 
the Eon between 2 and 3 Gyr BP, it appears that glacia- 
tions were altogether absent from the global surface. This 
is the crux of the faint sun paradox [Sagan & Mullen, 1972; 
Pollack, 1991; Sagan & Chyba, 1997]. 

A common solution to this apparent paradox is that larger 
amounts of the greenhouse gas (GHG) CO2 compensated for 
the cooler sun [Kuhn & Kasting, 1983; Kasting, 1993]. How- 
ever, some upper limits on the atmospheric partial pressure 
of CO2 (J3CO2) suggest that it was at most modest. For ex- 
ample, Rye et al. [1995] find pCO2^10 -1 ' 4 bar between 2.2 
and 2.7 Gyr BP, with indications that it could not have been 
much higher earlier [Sleep & Zahnle, 2001]. At these levels, 
the GHG warming could at most compensate for about half 
the reduction in the solar flux . Thus, it appears that CO2 
could have been important for warming the early Earth, 
but not sufficient to resolve the faint sun paradox by itself. 
Note however that theoretical considerations do suggest that 
PCO2 was higher in the early terrestrial atmosphere [Kast- 
ing, 1987] , such that it is likely that it did contribute at least 
some GHG warming. 

Because of these problems, it was suggested that other 
GHGs could have been important for warming the young 
Earth. In particular, it was suggested that small amounts of 
NH3 could have supplied the required GHG warming [Sagan 
& Mullen, 1972; Sagan & Chyba, 1997]. Although not im- 
possible, it is not easy to keep NH3 from irreversibly pho- 
tolyzing into H2 and N2. 

Another suggestion was that CH4 played the major GHG 
warmer [Pavlov et al., 2000]. This solution requires a long 
residency time of methane in the atmosphere, and prob- 
ably dominance of methanogenic bacteria. Unfortunately, 
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there are currently no reasonable bounds on the amounts of 
methane or ammonia in the early atmosphere. Thus, this 
type of solutions can neither be ruled out, nor proven at his 
point. 

Other effects are more certain but probably not large 
enough to solve the faint sun paradox. In particular, the 
young Earth rotated faster than today and it may have 
had less land mass. Employing a typical general circulation 
model (GCM), it was found that removing the land mass 
altogether could increase the temperature by 4°K, while a 
14 hr day would further increase the temperature by 1.5°K 
[Jenkins, 1993]. 

Another solution to the faint sun paradox involves strong 
solar winds, such that the earlier solar mass was higher. A 
10% more massive sun would have had a larger luminosity, 
enough to more than compensate for the early faint sun and 
also explain the observed Lithium and Beryllium depletion 
[Graedel et al., 1991]. However, solutions with a relatively 
gradual loss of mass from the sun, with a typical time scale 
of ~Gyr (as opposed to a mass loss time scale of a few 100 
Myr) are inconsistent with helioseismology, rendering this 
solution less probable [Guzik & Cox, 1995]. Moreover, a so- 
lar wind history can be reasonably reconstructed if the sun 
is compared to similar nearby stars at various ages. This 
gives a benign mass loss which is insufficient to resolve the 
paradox at any time of the terrestrial history [IToorf et al., 
2002]. 

Nevertheless, even a modest mass loss could potentially 
have interesting implications. The stronger wind was more 
effective at lowering the galactic cosmic ray flux (CRF). 
This, as we shall soon see, can present yet another possi- 
bility for helping resolve the faint sun paradox. It should 
however be stressed that it is not unlikely that the true so- 
lution is a combination of several of the above factors. 

Various empirical evidence [Tinsley & Been, 1991; Svens- 
mark, 1998; Egorova et al, 2000; Marsh & Svensmark, 2000; 
Todd & Kniveton, 2001; Shaviv, 2002a,b; Shaviv & Veizer, 
2003], theoretical analyses [Ney, 1959; Dickinson, 1975; Har- 
rison & Aplin, 2001; Yu, 2002], and experimental results 
[Harrison & Aplin, 2001; Eichkorn et al., 2002] link cosmic 
rays (CRs) to climate on Earth. The apparent CRF/climate 
link could originate through an effect that atmospheric ion- 
ization, which is governed by the CRF, may have on cloud 
formation. The empirical evidence is based on several sig- 
nificant correlations between the CRF and climate on Earth 
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over times scales ranging from days to geological time scales. 
The theoretical and experimental evidence describe vari- 
ous elements of an emerging physical picture relating atmo- 
spheric ionization, through the formation of condensation 
nuclei by charged ion clusters, to the formation of clouds. 

Even with the above evidence, this CRF/climate link is 
still a highly debated topic. This is primarily because a 
robust physical picture is still missing. For an objective re- 
view, the reader is encouraged to read Carslaw et al. [2002] . 
It should be emphasized that the analysis which follows rests 
on the assumption that this link indeed exists — an issue 
which we will address again in list of caveats in §5. Irre- 
spective of the outcome of this debate, since there is now 
enough evidence to seriously consider the possibility of CRs 
as a variable affecting the climate, we should consider the 
various implications of this link while bearing in mind that 
the underlying basis deserves extensive study. 

We begin in §2 with an estimate of the total warming ef- 
fect that the strong early solar wind has had, and compare it 
to the cooling induced by the fainter sun. In §3, we use the 
actual solar wind evolution (as deduced from nearby stars) 
combined with the variable star formation rate (SFR) in the 
Milky Way (MW) to calculate the evolution of the terrestrial 
temperature. We summarize our results in §4 and discuss 
them in §5. 

2. Estimating the effect of a stronger solar 
wind 

We assume henceforth that cosmic rays do have a cool- 
ing effect . Under this assumption, the stronger wind of the 
young sun had a warming effect. This is because the stronger 
wind necessarily lowered the galactic CRF reaching the in- 
ner solar system. A lower CRF would in turn translate into 
higher temperatures through the CRF/climate link. Thus, 
a stonger solar wind should have had a warming effect that 
acted to compensate the solar dimness. 

To quantify this effect, two basic numbers are required. 
First, it is necessary to know the ratio between changes in 
the CRF and changes in the global temperature. This will 
allow us to translate solar wind variations into temperature 
change on Earth. Second, if we wish to compare this ef- 
fect with the lower temperature associated with the faint 
sun, and also with possible GHG warming, we need to know 
how changes in the radiation budget translate into a global 
temperature change. 

The relation between radiative forcing and tempera- 
ture change is still an open and highly discussed question 
in climatology. For a blackbody Earth, this number is 
0.30°K/(W/m 2 ). Earth however does not behave as a black 
body, and various positive and negative feedbacks act to ei- 
ther amplify or decrease this number by a factor A. Values 
obtained in GCMs [IPCC, 1995] range between A = 2 and 
4.5. 

There are however debated claims that various negative 
feedbacks should be considered, yielding a much smaller A 
of 0.5 - 1.3 [Lindzen et al, 2001], or even ~ 0.2 [Ou, 2001]. 
This number was also obtained from analyses of actual tem- 
perature variations. For example, from the rather quick re- 
turn of the global temperature to its "average" after a series 
of volcanic erruptions, it was deduced that A ~ 0.5 with 
AJ>1 being inconsistent with the effects of volcanic driving 
[Lindzen & Giannitsis, 1998]. On the other hand, Gregory 
et al. [2002] find that a lower limit of A = 1.3 (at 95% con- 
fidence) can be placed by considering the increase of ocean 
heat content. Another constraint was obtained by Covey 
et al. [1996] who compared J3CO2 changes to temperature 
variations over geological time scales (in particular, the last 
glacial maximum, the Crataceous and early Eocene). The 



authors of this study found that A is consistent with the 
range of values obtained in GCMs. However, the latter two 
constraints were obtained while neglecting the possible forc- 
ing of the CRF and the nonthermal solar activity. Thus, 
if CRF/climate link is real, these two constraints are not 
applicable anymore. 

Once CRF variations are considered over geological time 
scales, a markable correlation is obtained between CRF as 
a driver and temperature change. This correlation, com- 
bined with the lack of correlation with pC02 variations, 
yields A ~ 0.75 (or A < 1.2 at 90% confidence and A < 2.2 
at 99% confidence) [Shaviv & Veizer, 2003]. Thus, if the 
CRF/climate connection is valid, which we assume here to 
be the case, we should consider A ~ 0.75 as the most prob- 
able value, but also A = 2.2 as an extreme upper limit. 

The same combined geological and astronomical analy- 
sis [Shaviv & Veizer, 2003] could also measure the relation 
between temperature change and relative CRF variations, 
which was parameterized through the form 

ATcrf « D [1 - (S /£ )«] , (1) 

where £ is the energy flux reaching the troposphere (to which 
the ionization is proportional), and £q is the flux reaching 
today. The power q relates changes in the ionization rate 
to changes in global temperature. Since the expected route 
is most likely through the effect that atmospheric ionization 
has on the formation of condensation nuclei, we should have 
q i [Yu, 2002; Harrison & Aplin, 2001]. It was also found 
that D = 10 ± 5°C, where most of the uncertainty is due 
to the inaccurate determination of the CRF variations from 
meteorites. We reduced here the upper range. Although a 
larger D is consistent with the meteoritic data, it is physi- 
cally inconsistent as it requires too large an albedo variation 
of Earth to be physically reasonable, as will soon be evident. 

When the solar system was young, the solar wind is ex- 
pected to have been significantly stronger. This implies 
that most of the CRF responsible for tropospheric ioniza- 
tion was blocked (i.e., the energy flux satisfied £ <C £o). 
Today, for comparison, the flux variations are only O(10%). 
Thus, we should expect an average global surface tempera- 
ture increase of AT w i nl i ~ D ~ 10 ± 5°K, when the wind 
was significantly stronger than today. 

Next, the decrease in solar luminosity as we go back to 
the end of the heavy bombardment era (at about 4 Gyr BP), 
is about ~25%. This is an inevitable result of all standard 
models 3 . The only astrophysical remedy is having a heavier 
younger sun. However, as mentioned in the introduction, a 
strong and slowly decaying wind that would allow the young 
sun to be bright for a long enough duration, is inconsistent 
with the winds of nearby stars [ Wood et al. , 2002] . 

The 25% reduction in solar flux corresponds to a decrease 
of ~60 W m~ 2 relative to the 240 W m~ 2 reaching the sur- 
face today. For A = 0.5, 0.75 or 2, the average global surface 
temperature decrease ATq associated with the fainter sun 
is ATq = -9°K, -13.5°K or -36°K respectively. Thus, for 
climate sensitivities consistent with the CRF/climate link 
(A<0.75), we find that the solar dimness can in principle be 
completely compensated for. 

Nevertheless, to actually obtain AT win d > |ATq|, wc 
require values of A and D which are consistent with ob- 
servations but are theoretically uncomfortable. For exam- 
ple, D = 15°K with A = 0.5 implies that CRs change the 
albedo by 28%, which is nearly impossible. The more rea- 
sonable values of D = 10°K and A = 0.75 yield a more 
realistic albedo change of 12.5%. In this case however, 
AT wind + ATq « — 3.5°K relative to today. This implies 
that the average Hadean and Archean temperatures would 
have been well above freezing. However, a full resolution 
of the faint sun paradox requires the temperatures at these 
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epochs to have been higher than today, since current day 
Earth has glaciations, while glaciations where absent from 
the young Earth (e.g., Crowell [1999]). 

With A = 0.75, the missing 3.5° K can be compensated 
with 16 W oT 2 . This could arise from a "modest" GHG 
warming. For example, 0.02 bar of CO2, which is about half 
of the present upper limits at ~ 3 Gyr BP [Rye et al., 1995; 
Sleep & Zahnle, 2001] would provide ~ 32 W m~ 2 of GHG 
heating 4 , ensuring that the early climate was warm enough 
to inhibit the appearance of glaciations. 



3. Solar wind evolution 
Temperature change 



To estimate the secular change in the solar wind modu- 
lation of the CRF, we first require the evolution of the solar 
wind parameters (i.e., mass loss rate and velocity). We will 
then translate this evolution to variations in the CRF reach- 
ing Earth. 

For the evolution of the mass loss rate M, we use the 
observational results of Wood et al. [2002] who found that 
the mass loss rate of nearby solar-like stars (main sequence 
stars of spectral type G and K) evolves as: M oc t M with 
H — —2.00 ± 0.52. Next, we take the velocity of the wind to 
be constant. This is because the terminal wind velocity de- 
pends primarily on the escape velocity from the sonic radius 
of the wind, which is expected to change only logarithmically 
with the mass loss rate (e.g., Lamers & Cassinelli [1999]). 
This implies that to first approximation we can take 8^0 
in a relation of the form v wind oc M s , but we leave S as a 
free parameter for generality. 

The energy loss suffered by CRs as they journey from 
the heliopause to Earth is through multiple scattering off 
the magnetic field of the solar wind. The average energy 
lost, for CRs with energies significantly larger than their 
rest mass, can be approximated with AE w Rv w i n d/(3A) 
[Perko, 1987], where R is the radius of the heliopause and 
A is a constant proportional to k, the diffusion coefficient 
of CRs in the solar wind. This expression is valid only if 
AE < £ > m c 2 and R > 1AU. Thus, for the calcula- 
tion of the energy loss we further require the scaling of both 
R and k with M. 

The value of R can be obtained by equilibrating the solar 
wind ram pressure at its termination Mv w i n d/R with the 
average pressure of the interstellar medium (ISM), which 
we ass ume to remain constant on average. This yields 

fiocVAf"— «M (1+S)/2 . 

The diffusion coefficient k depends on the effective mean 
free path £ arising from the Larmour gyration of the CRs 
in the disordered magnetic field of the solar wind. Since 
I oc B^ 1 , we have k ~ £c oc B^ 1 . Thus, the last scaling we 
seek is in the form B oc M 1 . It would allow us to relate M 
to changes in AE. 

At the radius of Earth's orbit, the magnetic field en- 
ergy density in the solar wind is typically 50 times smaller 
than the kinetic energy density of the wind. However, mag- 
netic field fluctuations should decay as r~ 3//2 [Jokipii,1973], 
(whether frozen-in or Alfvenic wave fluctuations), implying 
that at a few solar radii, where the solar wind is accelerated, 
the magnetic field energy density (oc r~ 3 ) is comparable to 
the kinetic energy density (oc r~ ). Although the current 
mass loss is probably not unique, we find equipartition be- 
tween the two energy densities at the acceleration region of 
the wind. Thus, it is reasonable to assume that this ten- 
dency for equipartition remains also for different M's. In 
such a case, B 2 r A ^ 2 oc pv 2 wini oc Mv win d at any given radius. 
This implies that 7 = (1 + 5)/ 2 if equipartition remains. 

The combination of these scalings yields AE m 
Afi/2+7+3,5/2 oc i A , with p, EE //(1/2 + 7 + 3(5/2). With 
nominal mass loss evolution (// w 2), constant wind veloc- 
ity (8 « 0) and aforementioned equipartition (7 « 1/2), we 
have fisi/iRi —2. 



The CRF spectrum (i.e., the differential number flux) in 
the ISM can be written as d^/dEisM = /ism = CE~ (p+1) 
for energies greater than about 1 GeV, where 3> is the CR 
number flux, C is a normalization constant, and p « 1.7. 
Furthermore, the differential number flux reaching Earth 
is related to fisM through f®(E) = fisM{E + AE). Since 
only CRs with an energy larger than the cutoff energy 
E c ~ 12 GeV can actually reach the troposphere, we find 
that the total energy £ > e c deposited in the troposphere is 



and predicted £ >Ec = 



f®E®dE® 



(2) 



C(£ e + AE)~ (P+1) E 9 dE t 



C(E c + AE)- p [^ T + — 
\p-l p 



More interesting is the energy flux at a time t (after the 
formation of the solar system) relative to the flux today, at 
to- We define e ee AE /E c . Since AE w 1.0 GeV [Perko, 
1987], we have e « 0.08, and also 



£>eM 

£>E C (to) 



(t/t ) 



P + e 



l + (t/U)-»e 
1 + e 



• (3) 



Eq. 1 for the temperature change induced by the wind 
evolution now becomes: 



AT 



. £ >Ec (t) S(t) Y 

£>E c (t ) S(t ) J 



D, 



(4) 



where S is the nearby SFR to which the CRF is 
proportional. For the SFR, we take the results of Rocha 
Pinto et al. [2000] who derived the MW SFR by analyz- 
ing the chromospheric ages of hundreds of late type dwarfs. 
The results of this comprehensive study are consistent with 
previous analyses, such as Barry [1988] who too used chro- 
mospheric ages, or Scalo [1987] who studied the mass dis- 
tribution of nearby stars. In particular, it appears that the 
SFR was higher than average between 2 and 3 Gyr before 
present, while it was below average between 1 and 2 Gyr 
BP. 

Interestingly, the MW SFR activity is correlated with 
SFR activity in the Large Magellanic Cloud (LMC). It ap- 
pears that there was a significant increase in the SFR in the 
LMC around 2 Gyr BP [Gallagher et al., 1996; Vallenari et 
al, 1996; Dopita et al, 1997] and that between 0.7 to 2 Gyr 
BP, the SFR in the LMC was below its average [Wester- 
lund,1990]. This correlation could be coincidental, but may 
reflect the tidal interactions between the two galaxies, dur- 
ing their perigalactica sometime around 0.2-0.5 Gyr BP and 
around 1.6-2.6 Gyr BP [Gardiner et al., 1994; Lin et al., 
1995]. Irrespective of the origin of the SFR variability, it 
modulated the CRF reaching the solar system, and presum- 
ably also climate on Earth. 

The fact that the moon was closer to Earth, and both 
were rotating faster [Zharkov, 2000] and that Earth proba- 
bly had a smaller land mass are not large effects by them- 
selves [Jenkins, 1993]. We will not try to model them accu- 
rately. Nevertheless, Jenkins [1993] has shown that no land 
mass and a 14 hr day will correspond to an average global 
temperature increase of typically 5.5°K. This however was 
obtained with GCM models which typically have A ~ 2.5. 
If, on the other hand, negative feedbacks exist to reduce A, 
the increase in temperature will be smaller accordingly. For 
A w 0.75, we should expect that the whole effect of no land 
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mass and a faster rotation would contribute only a ~ 1.6°K 
warming that gradually decreases with time between the 
formation of Earth and today. 

To mimic these two effects, we approximate the tempera- 
ture correction using AT corr (t) = A [(to — t)/to] 2.2°K where 
to = 4.5 Gyr is the time today. This is not a large correction 
at all, but unlike effects such as GHG warming by NH3 or 
CH4, which can potentially be much larger, the effects of 
land mass and rotation can be considered somewhat more 
quantitatively. 



4. Results 

Using the above physical components, we proceed to cal- 
culate the evolution of the average terrestrial temperature. 
The results are depicted in fig. 1 for various models which 
either include or exclude the solar wind/CRF effect using 
nominal values for the parameters, the effects of radiative 
forcing, and also the possibility of GHG warming by modest 
levels of pC02. 




Figure 1. Predicted temperature as a function of time before present for various models. Model (A) 
includes a radiative feedback of A = 2.2 and the GHG warming by modest levels of CO2 (by 0.01 bar of 
CO2 before 3 Gyr and exponentially decreasing afterwards to current levels). Model (B) includes no CO2 
contribution but A = 0.75. Model (C) has both CO2 and A = 0.75, but still has no CRF contribution. 
Model (D) includes CRF contribution using nominal values (fi = 2, E c = 12 GeV, AE = 1.0 GeV, 
D — 10° C). The short dashed line assumes the CRF reaching the outskirts of the solar system is constant 
and equal to the current day average. The solid line (no.l) assumes that the CRF is proportional to 
the SFR in the Milky Way [RochaPinto et al., 2000]. The dotted lines (no. 2) depict the temperature 
uncertainly associated with the SFR uncertainly error. The temperature variability associated with "short 
term" variability in the CRF, originating from the spiral structure [Shaviv, 2002a,b], is marked with the 
long dashed lines (no.3). These denote the temperature variability range. The actual temperature 
variability in the past 500 Myr [Shaviv & Veizer, 2003], which is associated with the variable CRF (as 
extracted from Iron meteorites [Shaviv, 2002b]), is marked with the short-dash-dotted line (no. 4). 



The figure recapitulates the main conclusion of 
§2. Namely, coupling the CRF/climate effect to the 
stronger wind of the young sun, results with a higher 
temperature. The effect is large, but if the nominal 
values are used it compensates only ~ |'s of the tem- 
perature reduction associated with the dim sun. Thus, 



either the CRF effect is larger than inferred from the 
nominal parameters, or more likely, there is some GHG 
warming. 

By estimating the actual evolution of the tempera- 
ture, we find additional conclusions to those obtained 
in §2. First, the effects of SFR variability, spiral arm 
passages and GHG warming are all comparable in am- 
plitude. Second, because of the different time scales 
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associated with the slowly varying SFR and the "fast" 
oscillations associated with spiral arm passages (over 
~ 10 8 yr), we find that at any epoch with a given SFR, 
a range of temperatures should exist from spiral arm 
variability. In other words, in addition to the statistical 
error on the average SFR [RochaPinto et at, 2000], we 
should consider that the location within the MW gives 
rise to additional variability 6 . Thus, not only is the av- 
erage important, but also the expected lower limit of 
the temperature range if we wish to compare with geo- 
logical evidence for the presence of glaciations. On the 
other hand, if we had a good record for the presence 
of extremely high temperatures, then the expected up- 
per limit of the temperature range would have been the 
important variable. 

If we study fig. 1 and consider models with the nomi- 
nal parameters of A w 0.75 and D w 10°K, we find that 
only during the past Eon (Neo-Proterozoic-Phanerozoic 
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Figure 2. A parameter study. Each panel describes 
models in which one parameter is changed from the nomi- 
nal model of fig. 1. The average global temperature is de- 
picted here with either a solid or dashed line (equivalent 
to the solid line of fig. 1). The expected range of vari- 
ability arising from spiral arm passages is depicted here 
with hatched shading (equivalent to the two-long dashed 
lines of fig. 1). Panel (A) describes the effect of pCO-2- 
The solid line is the nominal model with pC02 kept at to- 
days value while the dashed line has a saturation pC02 of 
0.02 bar before 3 Gyr BP. Panel (B) describes the effect 
of changing the mass loss evolution power law /i (within 
observed error range). Panel (G) describes the effect of 
the power q which relates the atmospheric ionization rate 
change to temperature change. Panel (D) describes the 
effect of changing the ratio between the average SFR in 
the past few 100 Myr and the current rate (nominal value 
is S = 0.7, Lavielle [1999]). Panel (E) describes the ef- 
fect of a changed climatic sensitivity to radiation forcing. 
Panel (F) describes the effect of a changed climatic sensi- 
tivity to CRF changes. The normalization is such that for 
today's flux, which is about 30% higher than the average 
over the past eon, the temperature is 15°C. 



period), and the period of 2-3 Gyr BP (late-Archean- 
Huronian) could there have been climatic periods dur- 
ing which it was cold enough as to produce tempera- 
tures lower than today and explain glaciations. This is 
consistent with the actual geological evidence, which 
indicates that only during the Neo-Proterozoic and 
Phanerozoic (past Eon), during the Huronean (2.2-2.4 
Gyr BP) and late Archean (2.9-3.0 Gyr BP) were glacia- 
tions present on Earth (e.g., Crowell [1999]). 

The exception is the first eon, which for the nomi- 
nal parameters should have been cold enough to have 
had glaciations (though not as cold as 2.5 Gyr BP). 
Since there is no evidence to support their existence, 
we return back to the point that even a modest GHG 
contribution, for example, by a modest 0.01 bar CO2 at- 
mosphere, can raise the temperature of the very young 
Earth such that with the reduced CRF, the temperature 
is kept above today's values. 

Fig. 2 is a parameter study showing the effects of 
changing various model parameters from their nominal 
values. It is apparent from the figure that there are two 
types of parameters. The first type are those which af- 
fect the total CRF/solar wind effect when the solar wind 
is completely switched off, or affect the way in which 
the effect of the dim sun translates into a temperature 
change. These parameters determine the gross temper- 
ature behavior. For example, A, D and maximum level 
of pC02 determine the temperature on the very young 
Earth. The second type of parameters are secondary in 
the sense that they only affect the evolution of the tem- 
perature between the young Earth and today. A good 
example is the power \x determining the temporal evo- 
lution of M. It does not affect the temperature in the 
extremely young atmosphere nor the temperature to- 
day, since the former is obtained in the limit of a very 
strong wind (irrespective of how strong it actually is), 
and the latter is simply given by the current M. 

Among the secondary parameters, the most impor- 
tant is the power q which relates the ionization rate to 
temperature fluctuations. The value q = i is theoreti- 
cally preferred [Harrison & Aplin, 2001; Yu, 2002]. A 
higher q would imply that CRF variability arising from 
a variable SFR and passages through spiral arms be- 
comes more prominent. Nevertheless, all the secondary 
parameters are not particularly important. 

Although future variability of the MW SFR cannot 
be predicted, we can predict the future evolution of the 
solar luminosity and the decrease in strength of the so- 
lar wind. Since the solar wind currently filters out only 
about 10% of the cosmic rays responsible for tropo- 
spheric ionization (with energies >10 GeV), the cool- 
ing effect of the solar- wind/CRF coupling is almost at 
its fullest. This implies that from now on, in the com- 
ing few billion years, the global temperature is expected 
only to increase on average. This is exemplified in the 
fig. 3. This is not to say that bursts of star forma- 
tion could reduce the temperature "temporarily", for 
example, when the LMC will be swallowed by the MW. 
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Time BP [Gyr] 

Figure 3. Past and future long term global temperature changes. The graphs depicted in this figure are 
the same as fig. 1, with a few exceptions, (a) "Short term" (~ 10 8 yr) temperature variability is not plot- 
ted, (b) Models with CO2 have pC02 = 0.04 atm before 3 Gyr BP (consistent with upper limit by Rye 
et al. [1995]), then decrease exponentially until today and remain at today's levels in the future, (c) The 
future temperature trend is plotted. Since the future SFR is unknown, the upper and lower temperatures 
associated with the actual SFR and spiral structure cannot be extended into the future. Instead, only 
the total range associated with both past SFR and spiral variability can be plotted. In the first ~ 5 Gyr 
of Earth's existence, the temperature increase due the increased luminosity is compensated for by the 
weaker solar wind and ensuing higher CRF. Since almost all the CRF (causing tropospheric ionization) 
currently reaches the inner solar system, this compensation cannot continue and the temperature will 
start to significantly increase. 



We expect that within about 1.5 Gyr, the long 
term average global temperature will increase to val- 
ues higher than anytime previously (i.e., >20°C), at 
which point the negative temperature excursions aris- 
ing from spiral arm passages and a higher SFR will not 
cool Earth down to current day temperatures, and no 
more ice-age epochs will ever appear. From about 2 Gyr 
into the future, even spiral arm passages and a higher 
SFR will not be able to reduce the global temperature 
to below the warmest temperature experienced in the 
past 4 Gyr. 

5. Summary 

We combined here the proposition that CRF varia- 
tions affect the global climate with the prediction that 
the early solar wind should have been stronger than 
today. This was shown to result with the following con- 
clusions: 

1. The stronger wind of the young sun was more ef- 
fective at blocking the CRF reaching Earth such that 
the CRF was significantly reduced when compared with 
the current flux. 



2. When coupled to the suggested CRF/temperature 
link, the lower CRF reaching Earth should have trans- 
lated into a higher global temperature. 

3. Using the apparent correlation between CRF vari- 
ability and global temperature over the past 550 Myr 
[Shaviv & Veizer, 2003], it is possible to quantify this 
effect. It was found that the associated warming was 
enough to significantly compensate for the faint sun and 
explain about ^ to |'s of the temperature increase re- 
quired to warm the young Earth to above the average 
temperature today. The remaining ^ to ^ could come 
from one of the previously suggested factors, such as 
GHG warming by modest levels of CO2 (pC02 of or- 
der 0.01 bar), which are consistent with the geological 
constraints. 

4. Besides explaining the secular climate trend, the 
evolution of the solar wind also reproduces, essentially 
without any free parameters, long term temperature 
trends on Earth once the variable SFR in the MW is 
considered. In particular, the CRF evolution explains 
why Earth has had glaciations in the past Eon and be- 
tween 2 and 3 Gyr BP. 

5. Although the "fine" details of a variable SFR and 
the passage through spiral arms cannot be predicted to 
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more than several 100 Myr into the future, it is possi- 
ble to extrapolate the decrease in solar wind and the 
increase in solar flux. It is found that from about 1.5 
billion years into the future, Earth is expected to be 
warmer on average than anytime in its previous history. 
This is because the CRF reaching Earth is currently 
about as large as it can be. 

These conclusions, however, do not come without 
caveats that are at least as important as the conclusions 
themselves. 

1. First and foremost, the aforementioned conclu- 
sions rest on the assumption that the CRF/climate con- 
nection is indeed real. Since this link still lacks a robust 
physical picture to support it (and is therefore still a 
subject of debate), it clearly deserves extensive study, 
especially in view of the repercussions that it has. For 
example, it bears consequences on the value of the ra- 
diative forcing, since only a low value is consistent with 
the CRF/climate link. 

2. Even if we were safe to assume that the CRF cli- 
mate connection is real, since the process is not fully 
understood, it is not trivial to guaranty that it also op- 
erated before the past Eon, during which periodic ice- 
age epochs existed. In particular, the CRF/climate link 
could depend on a biogenic agent which may have been 
absent during the Archean or early Proterozoic Earth. 
The link could also depend on atmospheric composition 
which undergone a few notable changes over the Eons. 

3. We have neglected in the present analysis ad- 
ditional effects that may be important but hard to 
quantify. These include for example effects of atmo- 
spheric composition, land mass characteristics (area 
and albedo), and even galactic migration of the solar 
system. 
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Notes 

1. This assumes that a change of lW/m 2 in the global radiation 
budget corresponds to the same global temperature change ir- 
respective of whether it is due to extra solar flux in optical and 
near-IR or the blocking of far-IR due to a higher level of pCC>2- 
This assumption is reasonable but not necessary. 

2. This can be seen in two ways. First, radiation budget mea- 
surements indicate that low altitude clouds have a net cooling 
effect [IPCC, 1995]. Thus, the positive CRF/low altitude cloud 
cover correlation [Svensmark, 1998) implies that a higher CRF 
will have a net cooling effect under the suspected CRF/climate 
link. Second, the long term positive correlation between CRF 
variations and ice-age epochs suggests that there is a nega- 
tive correlation between CRF and temperature [Shaviv, 2002a, 
Shaviv & Veizer, 2003]. 

3. It originates from the fact that as nuclear burning proceeds, 
the molecular weight of the solar core becomes progressively 
larger. 



4. We assume that the radiative forcing of CO2 on the atmosphere 
is F C 02 = 4W/m 2 [1.236 ln(c + 0.005c 2 ) - (ln(c + 0.005c 2 )] 
where co = 280ppm is the pre-industrial level of atmospheric 
CO2 [Hansen et al, 1998] 

5. Almost all supernovae in the MW are of types other than la. 
These originate from the death of massive stars which live less 
than ~ 30 Myr. Because this time scale is practically instan- 
taneous (relative to 4 Gyrs) and because CRs (at 10 GeV) 
originate from supernova explosions, the CRF will be propor- 
tional to the SFR (see references in Shaviv [2002b]). 

6. In the past Eon, for example, the average CRF was about 
~ 70% of the current CRF [Lavielle et al., 1999; Shaviv, 
2002b] but periods during which the solar system spent near 
spiral arms, the solar system was bathed with a flux higher 
than today, which presumably reduced the global temperatures 
enough to allow glaciations. 
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